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We report the forbidden 65^^5*0 — > 6s5d''Di magnetic-dipole transition amplitude computed 
using multi-reference Fock-space coupled-cluster theory. Our computed transition matrix element 
(1.34 X 10~*/ifl) is in excellent agreement with the experimental value (1.33 x 10"* Ms). This value 
in combination with other known quantities will be helpful to determine the parity non-conserving 
amplitude for the 6s^ ^So — > 6s5d^Di transition in atomic Yb. To our knowledge our calculation 
is the most accurate to date and can be very important in the search of physics beyond the standard 
model. We further report the 6s6p^Po — > 6s6p^Pi and 6s5d^Di — > 6s6p^Po transition matrix 
elements which are also in good agreement with the earlier theoretical estimates. 
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I. INTRODUCTION 

The highly forbidden 6s^ ^5*0 — > 6s5d ^Di magnetic-dipole (Ml) transition amplitude in ytterbium (Yb), 
a key quantity for evaluating the feasibility of parity non-conservation (PNC), has recently been measured 
by Stalnaker et al. using Stark-interference experiment. The electric-dipole (El) matrix element for 
6s^ ^5*0 — > 6s5d '^Di transition in Yb is forbidden because of its s — d nature. The forbidden Ml transition 
amplitude mentioned above is therefore the key quantity to explore the feasibility of the PNC study for this 
transition in Yb. Accurate determination of the Ml transition amplitude, which is strongly suppressed in 
nature in the absence of external fields, can be used together with the large PNC- and moderately large 
Stark- induced El amplitudes to understand PNC studies in neutral Yb. Strong configuration mixing and 
spin-orbit interaction in both the upper and the lower states give rise to a non-zero 6s^ ^5*0 — > 6s5d ^Di 
transition amplitude [H, Q- Surprisingly, despite its tremendous importance in PNC experiments, only a 
rough theoretical estimate (|^(M1)| ^ lO"'*/!^) is available in the literature for this transition. PNC in 
atoms arises from the neutral weak interactions and are considerably enhanced in heavy atoms. Combining 
the high precision experiments and theoretical calculations of PNC observables, it is possible to extract the 
nuclear weak charge Any discrepancy of its value with the one obtained from the standard model (SM) 
of particle physics could possibly reveal the existence of new physics beyond the SM. 

The ground and excited states of closed shell ground state systems Hke Yb are, in general, multi- 
configurational in nature, and hence, an accurate description of these states requires a balanced treatment 
of non-dynamical or configuration mixing and dynamical electron correlation effects (this will be more clear 
by studying the energy levels in figure [T]). It is, therefore, imperative that these systems must be treated 
with methods which are combinations of the configuration interaction (CI) and many-body perturbation 
theory (MBPT), such as multi-reference many-body perturbation theories (MR-MBP_T) J, |5, 6, I, 8, 9], 

[M, m M etc. 



multi- reference Fock-space coupled-cluster (MR-FSCC) theories and/or it variants [lOl. Illl. 
The state-of-the-art MR-FSCC is an all-order approach and is capable of providing reliable estimates of pre- 
dicted quantities. In this paper, we employ the MR-FSCC method to compute the magnetic-dipole transition 
amplitude for 6s^ ^5*0 — > 6s5d ^Di transition in Yb using four-component relativistic spinors. The resulting 
value of this magnetic-dipole transition matrix element in atomic Yb is 1.34 x 10"'* /x^, which differs by less 
than one percent from the experimental value. In addition, we have also calculated the 6s6p^Po — > 6s6p ^Pi 
Ml transition transition amplitude in Yb which plays crucial role in the measurement of the PNC induced 
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Figure 1: Energy levels of the ground and low lying excited states of Yb. The energies (in cm^^) are given with 
respect to the ground state and are obtained from the NIST database jlBl- Electric dipole (allowed) and magnetic 
dipole (forbidden) transitions are represented by 'blue' (solid) and 'red' (dashed) lines respectively. This diagram 
helps us to understand the requirement of a multi-reference theory to describe the atomic states of Yb. 

electric-dipole amplitude [H] . This is the first time any variant of coupled-cluster theory has been appHed 
to determine the Ml transition amplitude of Yb. A precise determination of this quantity ensures not only 
the power of the theory but also for the experimental uncertainties. To our knowledge no such theoretical 
results are available for magnetic dipole transitions in Yb. 

The structure of this paper is the following : section U] describes the physical relevance of the problem. 
Section [TTl provides a brief outline of the multi-reference Fock-space CC (MR-FSCC) theory for two-electron 
attachment processes that is used to compute the Ml transition elements between the ground ^5*0 and excited 
^Di state. Section Hill contains the results of our calculation with an in-depth discussion. Finally in section 
ITVl we conclude and highlight the findings of our paper. 

II. FOCK-SPACE MULTI-REFERENCE COUPLED-CLUSTER (MR-FSCC) THEORY FOR 
TWO-ELECTRON ATTACHMENT PROCESSES 

In MR-FSCC method [HI, [H, [13, H H llOl , the self-consistent field (SCF) solution of the Hartree- 
Fock (Dirac- Fock in relativistic regime) for the TV-electron closed shell ground state $hf/df is chosen as 
the vacuum (for labeling purpose only) to define holes and particles with respect to $hf/df- The multi- 
reference aspect is then introduced by subdividing the hole and particle orbitals into active and inactive 
categories, where different occupations of the active orbitals will define a multi-reference model space for 
our problem. We call a model space to be complete if it has all possible electron occupancies in the active 
orbitals, otherwise incomplete. The classification of orbitals into active and inactive groups is, in principle, 
arbitrary and is at our disposal. However, for the sake of computational convenience, we treat only a few 
hole and particle orbitals as active, namely those are close to the Fermi level. The classification of orbitals 
is depicted schematically in FiglSUa) . Diagrammatically, active holes and particles are depicted as solid lines 
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with double arrows and the corresponding inactive lines are designated by dotted lines with single arrow. The 
orbitals which can be both active and inactive are designated by solid Hnes with single arrow (see FiglUb)). 

We designate by \I>°''^''^ a model space of k-ho\e and Z-particle determinants, where in the present instance 
(Yb^'^ + 2e — > Yb), fc = and I ranges from to 2. Generally, any second quantized operator has fc-hole 
and Z-particle annihilation operators for the active holes and particles. For convenience, we indicate the 
"hole-particle valence rank" of an operator by a superscript (fc, I) on the operator. Thus, according to our 
notation, an operator ^C^'') will have exactly fc-hole and ^-particle annihilation operators. 



active 



inactive 
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Figure 2: (a) Schematic depiction of the classification of particle and hole orbitals into active and inactive categories, 
(b) Diagrammatic representation of hole (|), particles (j), active particles (double up arrow), inactive holes, and 
particles (dashed down/up arrow), (c) Diagrammatic representation of S'*-"'''-' (T), S''"'^', and S'"'^' cluster operators. 



We now describe the type of ansatz used to derive the MR-FSCC equations for direct energy difference 
calculations in two-electron attachment processes. The Hartree-Fock/Dirac-Fock function $hf/df is denoted 
by ^fC'O) and the inactive hole and particle orbitals (defined with respect to $hf /dp) are labeled by the indices 
a, 6, c, ■ • ■ and p,q,r,- ■ ■, respectively. The corresponding active holes and particles are labeled by the indices 
a,/3, 7--- and u,v,w---, respectively. Note that there will be no active holes (particles) for two electron 
attachment (detachment) processes. The cluster operator correlating the iV-electron ground/reference state 
is denoted in our notation by S'^^'*^-' which can be split into various rt-body components depending upon 
the various hole-particle excitation ranks. The cluster operator upto 2-body (first two diagrams of 

Figlljc)) can be written in second quantized notation as, 

^(0,0) ^ ^(0,0) ^ ^(0,0) ^ _ J2iP\^f-'^\-)H-^} + \ E {P<l\st'^\ab){alala,a^} + ■■■ (1) 

p,a a,b,p,q 

where (a) denotes creation (annihilation) operator with respect to 3'hf/df and {• ■ •} denotes normal 
ordering. It should be noted that S'(°'°) cannot destroy any holes or particles; acting of $hf/dF) it can only 
create them. 

For (A^+1) electron states the model space consists of zero active hole and one active particle (fc = 0, Z = 1) 
and hence according to our notation the valence sector for {N + 1) electron states can be written as (0,1) 
sector. We introduce an wave operator fl which generates all valid excitation from the model space function 
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for {N + 1) electron states. The wave operator fl for the (0,1) valence problem is given by 

f} = {exp(5("-0)+5(0-i))}. (2) 

In this case the additional cluster operator S^^'^^ must be able to destroy active particle present in the 
(0,1) valence space. Like 5"*^°'°^ the cluster operator S''-"'^^ can also be split into various n-body components 
depending upon hole-particle excitation ranks. The one- and two-body S'^^'^^ (3rd and 4th diagram of 
Figl2|c)) can be written in the second quantized notation as 

^(0,1) ^ ^(0.1) ^ ^(0,1) ^ . . _ ^(p|4o,i)|^^{^t,„} + 1 ^ {pq\sf^'^\ua){alala,a^} + ■■■ (3) 

where u denotes the active particle which is destroyed. 

Similarly, for {N + 2) electron states (two-electron attachment processes) the model space consists of zero 
active hole and two active particles (k — Q,! — 2) and the valence sector may be written as (0,2) sector. 
In this case, the additional cluster operator must be able to destroy two active particles and this may be 
designated by S^^''^\ The total wave operator for the (0,2) problem is then given by 

n = {exp(5("^o)-H5(°'i)+5(°'2))|^ 
A typical two-body S'j"'^'' operator (5th and 6th diagrams of FiglJ] (c)) may be written as 

"^2"'^^ = X! {P<l\s2^'^^\u'v)Wpa\a^au}, (5) 
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where u and v denote active particle which are destroyed. Note that orbitals p and q both cannot be active 
at the same time. We further emphasize that under two-body truncation scheme S''^°'^^=0, if all the particle 
orbitals are active. 

In general, for a (fc, I) valence problem, the cluster operator must be able to destroy any subset of k- active 
holes and I- active particles. Hence, the wave operator fl for (fc, I) valence sector may be written as 

f] = {exp(^(^-''))}, (6) 

where 

k I 

S{k,i) ^ J2 ^S'(™."). (7) 

m— n— 

To compute the ground to excited state transition energies and Ml transition element(s) of Yb, we begin 
with the Dirac-Coulomb Hamiltonian (H) for an A^-electron atom which can be written as 

AT N 2 

- 1 ■ ^ 

with all the standard notations often used. The normal ordered form of the above Hamiltonian, relative to 
the mean field energy, is given by 

H = H- = H-Edf^ E^^I^IJ') {"I'^^ + i ^ ■ (9) 

Here 

mki) = {i3\—\ki)~{ij\—\ik), (10) 

Eup is the Dirac-Fock energy and f is the one-electron Fock operator. 
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We define the exact wave function for (fc, I) valence sector as 

where 

i 

The functfons $f ''^ in Eq.^ are the determinants included in the model space are the 

corresponding coefficients. Substituting the above form of the wave-function (given in Eqs. ifTTj) and (fT2| l 
in the Schrodinger equation for a manifold of states = we get 

i/f](^^C,:|<i>f'))j ^E,n(j2C^\^l''^?j, (13) 

where is the «-th state energy. 

Following Lindgren [ll|, Mukherjee Lindgren and Mukherjee [Hi, Sinha ei al. [Hi and Pal ei al. [2q| . 
the Fock-space Bloch equation for the MR-FSCC may be written as 

where 

and pC^'') is the model space projection operator for the (fc, I) valence sector (defined by J^i c|'°''^$j-'^''^). For 
complete model space, the model space projector p(*^'') satisfies the intermediate normalization condition 

p{k,l)^p{k,l) ^ p(kl) ^ ^j^g-j 




(b) 

Figure 3: Diagrammatic representation of -ff^g^' (figure a) and H^g^' (figure b). The one- and two-body dressed 
operators of H are represented by circle and dashed lines, respectively. Exchange diagrams are not shown here for 
convenience. 

At this juncture, we single out the cluster amplitudes and call them T. The rest of the cluster 

amplitudes will henceforth be called S and are shown in Fig. [51 The normal ordered definition of enables 
us to rewrite Eq.|[7l) as 



n = exp(T){exp(S')} = exp(T)ri^ (17) 

where flv represents the wave-operator for the valence sector. 

To formulate the theory for direct energy differences, we pre-multiply Eq. (fT4|) by exp(— T) and get 

Ttn^pC^'i^ = VI^P''^'^^ H^^^'^ P^^'^^ , V(fc,Z) ^ (0,0) (18) 

where H ~ cxp(— T)iJ exp(T). Since H can be partitioned into a connected operator H and i?rcf/gr 
electron closed-shell reference or ground state energy), we likewise define Hcs as . 



(k,i) 

eff 



H, 



(k,i) 



eff 



i?g,, v(fc,0 7^(o,o). 



(19) 
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Figure 4: Diagrammatic representation of (^j ' ' ). The one-body dressed operator Ml = 

exp(T^)Afl exp(r) is represented by line with circle. Exchange diagrams are not shown here for convenience. 

Substituting Eq.ljlQp in Eq. lfTS]) we obtain the Fock-space Bloch equation for energy differences: 

i/f^^pC^-') = f7„p('='')if<^g.''V('^^'), V(fc, I) ^ (0, 0). (20) 

Eqs. lfT4|) and (|20l) are solved by the Bloch projection method for k = I = and fc = 0, ^ 7^ 0, respectively, 
involving the left projection of the equations with pC^'') and its orthogonal complement Q^'''^^ (p(fe.O -(- 
= to obtain the effective Hamiltonian and the cluster amplitudes, respectively. At this point, we 
recall that the cluster amplitudes in MR-FSCC are solved hierarchically through the subsystem embedding 
condition (SEC) (itI . [2iI | which is equivalent to the valence universality condition used by Lindgren [ll| in 
his formulation. For example, in the present application, we first solve the MR-FSCC for fc = Z = to obtain 
the cluster amplitudes T. The operator H and H^^'^^ are then constructed from this cluster amplitudes T 
to solve Eq. (|20|) for fc = 0, Z = 1 to determine 5'^*''^' amplitudes. The effective Hamiltonian for (0, 1) Fock 
space (represented diagrammatically in FigEl), constructed from H, T and 5"^°'^^ is then diagonalized within 
the model space to obtained the desired eigenvalues and eigenvectors. The diagonalization is followed from 
the eigenvalue equation 



jJ(o.i)^(o,i) ^ (21) 

where 

jj(o.i) ^ p(oa)[^+f5^]p(o,i). (22) 

The expression HS^'^'^'^ in Eq. (|22|) indicates that operators H and 5'*^°'^^) are connected by common orbital(s). 

The MR-FSCC equations for (0, 2) sector are then solved to determine where the cluster amplitudes 

from the lower valence sectors behave as "known quantities". The effective Hamiltonian for the (0, 2) Fock 
space constructed from H, T, S^°'^^ and 5''°'^^ is then diagonalized to get the desired roots by using the 
equation 

^(0,2)^(0,2) ^p(0,2)^_ (23) 

where 



^(0,2) ^p(0,2)[^^^5(04)+l'^5(0,l)5(0,l)^^5(0,2)]p(0,2)_ (24) 
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It is worth noting that the eigenvalue and eigenfunctions for the (0, 1) valence sector are by-products of 
MR-FSCC for the (0, 2) valence sector with no additional computation. Once the cluster amplitudes are 
known, the magnetic-dipole matrix element between the two states can be computed using the following 
expression 

(Finalstate(/)|Afl|Initialstate(i)) = ^ (25) 

' (0,2)|,.,(0,2).,,.,(0,2)|,.,(0,2) 



where and l^-f" ^) are the exact initial and final states, respectively. With aid of fi, the valence 

universal wave operator, Eq.lfTSll can be further simplified to 



($f '^1(1 + 5t)ertgT(i + 5)|$("'2)^^^(o,2)|^^ ^ 5t)erte^(l + S)\^f^^) 

where MT = exp(rt)Ml exp(r) and S = S^^'^^ + S^°^'^\ 

The single particle reduced matrix elements for the Ml transition is given by. 



{njW ml ||«,> = A (j-^ll c(i) ll^v) X 1 ji(fcr) (T'/Q. + Q/T'.) dr. (27) 

Here j's and k's are the total orbital angular momentum and the relativistic angular momentum quantum 
numbers respectively; k is defined as uja where oj is the single particle difference energy and a is the fine 
structure constant. The single particle orbitals are expressed in terms of the Dirac spinors with Vi and Qi 
as the large and small components for the zth spinor, respectively. The angular coefficients are the reduced 
matrix elements of the spherical tensor of rank m and are expressed as 



('^/llCf ) 11'^.) = (-l)-'>+^/y(2j/ + l)(2j- + l) 



3f 

I 

2 



TO 





(28) 



with 



n 1 \ ( 1 iilf + m + k even , , 

4Z/,m,/.) = | 'otherwise (29) 

and Vs being the orbital angular momentum quantum numbers. When kr is sufficiently small, the spherical 
Bessel function jn(fcr) is approximated as 

, (fcr)" (fcr)" , , 

J"(^")"(2;^=1.3.5....(2n+1)- ('"^ 



III. RESULTS AND DISCUSSIONS 



The magnetic (Ml) and electric-dipole (El) transition matrix elements of Yb are computed using 
37s33p28dl2/5g GTOs with ao = 0.00525 and /3 = 2.73 (geometrical basis with a; = ao/?'~^). [High 
lying unoccupied orbitals are not included (kept frozen) in CC calculations.] The reference space for ex- 
citation energy and associated properties is constructed by allocating 6s valence electrons of Yb among 
6s7s8s6p7p5d6d valence orbitals in all possible ways. The basis and reference space used in this calculation 
is exactly same as that employed in an earlier communication by one of the author for transition energies, 
ionization potential and hyperfine matrix element calculations. We have considered that the nucleus has a 
finite structure and is described by the two parameter Fermi nuclear distribution 



1 + exp((r — c)/a) 
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Table I: Theoretical and experimental magnetic dipole transition matrix elements (in Bohr magneton ^b) of Yb. 



Initial State Final State This work Expt. /Theory 

Gs" 'So 6s5d^Di 1.34 x lO"* 1.33 x lO^^fl] 
6s6p^Po GsGp^Po 0.12 0.13(16] 



where the parameter c is the half charge radius and a is related to skin thickness, defined as the interval 
of the nuclear thickness in which the nuclear charge density falls from near one to near zero. The en ergy 
levels of Yb and Yb^ are not reported here as those have already appeared in the previous work [2a |. 
The magnetic-dipole transition matrix elements in Yb computed using MR-FSCC method agree well with 
experiment and with other available theoretical calculations (see Table HI) The present result for |^(M1)| 
for Gs^^Sq — > 6s5d^Di transition differs by less than one percent (< 1%) from the experimental value. 
Our calculation further shows that the major contribution to |A(M1)| comes from StO'i) (5.(0.2) contribution 
is only 1%). At this juncture, we emphasize that the random phase approximation (RPA) and the second 
order multi-reference many-body perturbation theory (MR-MBPT(2)) estimate this quantity (|^(Afl)|) to 
be 0.68 X IO'^'/xb and 0.98 x 10~''/u_b respectively. These large deviations (~ 49% for RPA and ~ 25% for 
MBPT) in the perturbative estimate clearly demonstrates the importance of higher order correlation effects. 

In addition to the trasition matrix element (^6s5d^Di\ Ml\6s^ ^ So) , we also report the 6s6p^Po — > 
GsGp^Pi Ml transition amplitude in Yb, which plays an important role in the measurement of PNC induced 
electric-dipole amplitudes We briefly outHne its relevance as the details are available elsewhere 

The PNC-induced electric-dipole transition amplitude 74(i?l)pNc is given by 



A{El) 



PNC 



(6s6p ^Pi 
(5^3/2651/2 



ez |6s6p'^Po) 
|5d3/26pi/2> 



(32) 



where is the PNC weak interaction Hamiltonian in the non-relativistic limit, e is the electronic charge, b 
is a coefficient that describes the configuration mixing amplitude and angular mixing coefficient, and AE is 
the energy separation between the 6s5d^Di and 6s6p^Pi states The mixing coefficients of the 6s5d^Di 
and 6s6p^Pi states by the weak interaction is given in Ref.[3|. We have also determined the matrix element 
(6s5d^Di\ ez\6s6p^Po) which turns out to be 2.52 a. u. This value provides a step forward towards the 
determination of A{E1)pnc amplitude in Yb and in the search of physics beyond the standard model. 



IV. CONCLUSION 



We have computed the highly forbidden magnetic-dipole transition matrix elements for 6s^ ^Sq — > 
6s5d^Di and 6s6p^Po — > 6s6p^Pi transitions in Yb using the Fock-space multi-reference coupled-cluster 
(MR-FSCC) method. The values of the magnetic-dipole transition matrix elements presented here are the 
most accurate theoretical estimates to date and are in accord with the experimental value. We have also 
evaluated the (^6s5d^Di \ ez |6s6p^Po) matrix element, which can be combined with other known quantities 
to determine the PNC amplitude for the 6s^ — > Gsbd^Di transition in atomic Yb. To our knowledge 
this the first time any variant of coupled-cluster theory is applied to determine this quantity, which is ex- 
pected to be useful to experimentalists in this area and in the search of any new physics beyond the standard 
model. 
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